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Abstract 

The contribution of Lorentz non-scalar operators to finite temperature correlation func- 
tions is discussed. Using the local duality approach for the one-pion matrix element of a 
product of two vector currents, the temperature dependence of the average gluonic stress 
tensor is estimated in the chiral limit to be (E 2 + B 2 )y = ^bT A . At a normalization point 
fi = 0.5 GeV we obtain b « 1.1. Together with the known temperature dependence of 
the Lorentz scalar gluon condensate we are able to infer (E 2 )t and (B 2 )t separately in the 
low-temperature hadronic phase. 
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Correlators of currents with the quantum numbers of hadrons are known to be useful to 
obtain information about the masses and couplings of hadrons; they are employed in the 
QCD sum rule approach and in lattice calculations. In both approaches the correlators are 
considered at large Euclidean distances or imaginary times where the dominant contribution 
comes from the lowest state with the corresponding quantum numbers. QCD sum rules give 
predictions also for form factors and structure functions of hadrons. (For a recent review of 
applications of QCD correlation functions see ref. [|TJ.) 

In recent years there has been increasing interest in finite temperature QCD and hadronic 
physics due to the expectation that at high enough temperatures the QCD vacuum, spec- 
ified by nonperturbative condensates of quark and gluon fields, will "melt" and undergo a 
transition to a quark-gluon plasma. Melting is usually understood in the sense that chiral 
symmetry restoration and deconfinement take place. The former means that with increasing 
temperature quark condensates evaporate, while the latter means that hadrons do not rep- 
resent stable degrees of freedom. It was shown by Leutwyler and his collaborators |2| using 
the chiral Lagrangian approach that the quark condensate indeed decreases with rising tem- 
perature. From the usual QCD sum rules at T = it is well known that the properties of 
hadrons are, to a large extent, determined by nonperturbative quark and gluon condensates 
||. Naturally, a large number of papers were devoted to the generalization of QCD sum 
rules to finite temperature in attempts to relate the temperature dependence of the hadronic 
spectrum to the temperature dependence of the condensates (see, e.g. f|, |, ||). In this case 
the vacuum average of the product of currents becomes the Gibbs average over the thermal 
ensemble. To calculate the Gibbs average one must choose a basis for the states. As argued 
in refs. |5], [7] at temperatures which are much less than the energy scale of confinement the 
appropriate basis is that of hadronic states, rather than the quark-gluon basis used in early 
papers on the subject (see, e.g. ref. Q). Using this basis it was also shown H that at low 
T the thermal correlators are expressed as a mixture of zero-temperature correlators with 
different parity. It is also clear that if the operator product expansion (OPE) is applied to a 
thermal correlator then the temperature dependence appears only in the matrix elements of 
the operators (condensates), the coefficient functions being obtained through a perturbative 
calculation at T = 0. QCD sum rules at low temperature were recently reexamined along 
these lines in ref. [f|.[] At high temperatures, corresponding to the quark-gluon plasma, the 
calculation of thermal correlators should be performed in a basis consisting of quark and 
gluon states. In this case the perturbative temperature-dependent parts of the condensates 
due to quarks and gluons from the thermal ensemble may be included in the coefficient 

■"■We thank T. Hatsuda for drawing our attention to this paper. 
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functions ||. 

Thus the QCD sum rule method, understood as a tool to get information about the 
imaginary parts of correlators via analyticity, seems to be tractable both at very low and very 
high temperatures, but not in the region of a phase transition where a drastic rearrangement 
of the spectrum takes place. 

An additional feature of finite temperature sum rules is the appearence of new conden- 
sates due to Lorentz non-scalar operators; these were, of course, present in the OPE, but 
gave zero contribution when averaged over the vacuum. At finite temperatures Lorentz in- 
variance is broken and these operators should contribute |], |J. The same applies to the case 
of finite density ||. However, each of these new condensates is an unknown nonperturbative 
parameter. In principle they may be fixed from the physical spectral densities of the corre- 
lators, just as in the zero temperature case the now well-established condensates were fixed 
by the hadronic spectrum. 

Consider the correlator of two isovector vector currents at finite temperature T and 
euclidean momentum q, where T 2 <C Q 2 = —q 2 and Q 2 > 1 GeV 2 : 

i [ *xe*^(n\Tjp(x)j v (0)e^W\n) = {g,uq 2 -q,qu)CMT) + uy v C 2 {q,T) , (1) 

J n 

where T denotes a time-ordered product, = -(uj^u — d^^d), u* = — (u ■ q)q^/q 2 is the 
transverse part of the heat bath four-velocity and Vt = — T \og{J2n{ n \ e ~ H ^ T \ n ))- Eq. (0) is 
the most general expression compatible with conservation of the vector current. The Lorentz 
invariance breaking term proportional to utul, must be absent at T = 0. This means that 
C 2 (g, T) goes to zero as T — > 0, while Ci(q, T) becomes the usual zero temperature correlator. 
Notice that eq. (|I|) may be considered to be the amplitude for forward scattering of a virtual 
photon by the heat bath. Then the imaginary parts of C% and C 2 are the structure functions 
of deep inelastic scattering of leptons by the heat bath (v? is similar to the transverse 
component of the target momentum, p M — (p ■ q)q^/q 2 ). 

At low T, when the contributions from all particles except pions are exponentially sup- 
pressed in the Gibbs average, the functions C\ and C 2 may be estimated by expanding in 
the density of thermal pions. In the first order of this expansion only matrix elements over 
one-pion states are taken into account. This approximation was made in ref. [U for C±. 
The one-pion matrix elements were estimated via PCAC and current algebra. It was shown 
that C\ and its counterpart from the axial channel are given by T-dependent mixtures of 
their zero temperature values and, as a result, the corresponding screening lengths tend to 
converge with increasing temperature. 

The purpose of the present letter is to estimate C%. Let us start from the one-pion matrix 



2 



element in the chiral limit 

i I d^e^(n(p)\Tj,(x)UO)\rr(p)) , (2) 



where we assume j) ~ T < Q, since eq. (|j) is to be integrated over p with Bose occupa- 
tion probabilities. If AtlAt2 ... /Jn is an operator of Lorentz spin n, then the matrix element 
(n(p)\0 

ti 1 ti 2 ...tJ.n\ 7r {p)) P/j.iP/j.2 ■ ■ - Piini an d cannot be reduced via PCAC to a vacuum matrix 
element. It is clear that at low temperatures, T <C Q, the main contribution to Ci comes 
from operators of lowest spin, namely spin 2. In leading twist there are two spin 2 operators 
which are related to the energy-momentum tensor: 

Q li m = ¥^^ D ^<l + Tfn2 D tix ( i) , q = u,d } s... 

nG /~ia /~iaa 1 s~ta /~<aaf3 jn\ 

where _D M is the covariant derivative. Graphs which correspond to the contributions of these 
operators to the matrix element in eq. (|2|) are shown in fig. 1. If the normalization point 
for the operators is taken to be n 2 = Q 2 , then the operator #^ lAl2 does not contribute to the 
OPE in the leading log approximation, and the contribution of twist 2, spin 2 operators to 
eq. (0) involves 

^(*(p)K + *J,k(p)> = ^(p)K - o%Hp)) . (4) 

Here we neglected the contributions of heavy quarks. The matrix element of the total 
energy-momentum tensor is (7r(p)|0*^|7r(p)) = 2p il p v (the states are normalized such that 
(7r(p)|7r(p')) = (2ir) 3 2E5^ (p — p')), while the matrix element of the gluon energy-momentum 
tensor 

{*(p)\0%,W)) = bp ll p v (5) 

contains an unknown constant b. This constant is related to the matrix element of the energy 
density of the gluon field 
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v (7r(p)|E^ + B>(p)>, =Q . (6) 

Note that b depends on the normalization point, //, in the operator product expansion. This 
dependence will be discussed later. 

Let us try to estimate b within a quark-hadron duality approach, saturating the amplitude 
of eq. (D by hadrons, (7r|Tj M (x)j v (0)|7r) = J2n( 7! '\j^)\n)(n\j u (0)\iT). Focussing on spin 2 
contributions to C2, we then have 

2-6 , c 1 [°° p(s)F 2 (Q 2 ) 
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where F n (Q 2 ) is the part of the form factor {^{p)\jp\n{p + q)) proportional to p p and p(s) is 
the spectral density in the s-channel. The states \n) are normalized as in eq. (BJ), the n-state 
contribution to p(s) being 7iS(s — m 2 ). On the l.h.s. of eq. (|7|) the term c/Q 4 denotes the 
contribution of three different spin 2, twist 4 operators whose individual contributions 



cannot be separated. The constants b and c are considered as parameters to be fitted. The 
ellipsis in eq. (0) corresponds to spin 2 terms of higher twist. Note that eq. (|7|) is just 
the sum rule for the second moment of the deep inelastic structure function, Jq 1 ^(x, Q 2 )dx, 
divided by Q 2 . It is valid in the asymptotic region, Q 2 — > oo, with all higher states in the 
s-channel equally important in this region. Our goal here is to see whether eq. (|7]) can be 
satisfied in a region of intermediate Q 2 ~ 1 GeV 2 where the r.h.s. may be approximated by 
the contribution of a few low-lying statesQ. 

First consider the case of charged pions. The lowest states in the s-channel are the 7r and 
ai(1260) mesons. Assuming p-dominance for the form factors (which is known to be a good 
approximation for the pion form factor up to Q 2 ~ 2 GeV 2 ), (ir\j p \n) = — e^(iTp\n)(Q 2 + 
m 2 ) -1 , where is the p-meson polarization vector and g 2 /4:ii ~ 2.9. We obtain for the r.h.s. 
ofeq. © 



3m 4 



(Q 2 + m 2 ) 2 



Here we used 



Q 4m 2 ai g 2 p {Q 2 + mlj L 

(K + (p)p°(q)W + (p + q)) = 9 P ,n e p * ■ (2p + q) (9) 



and 

i(n + (p)p°(q)\at(p + q)) = g aifm s p * ■ e a ^ + h ai(m e p * ■ (p + q) e a ' ■ p . (10) 



The notation corresponds to that of ref. [0. Note that g piT1T = g p within the p-dominance 
approach. The couplings g aiP n and h aipn cannot, of course, be determined from the ai width 



alone. To this end we use an effective chiral Lagrangian with spin 1 mesons [12, 13, 14 



In this approach the constants in question are expressed in terms of parameters of this 



2 In ref. JlT[ the transverse photon structure function in the region of intermediate x was calculated 
starting from the WW four-point correlation function, using the OPE in the photon virtuality p 2 and 
extrapolating to p 2 = 0. One could think of doing the same thing for the pion structure function, starting 
from the AW A correlator. It can be shown, however, that just as in the case of the longitudinal photon 
structure function, there are difficulties in the extrapolation to on-shell pions. The AW A box diagram also 
cannot be used, via a triple dispersion relation, to model the continuum contribution to the real part of the 
forward scattering amplitude in the usual manner because of the zero momentum transfer in the t-channcl. 
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Lagrangian which are fitted to reproduce masses and widths. The Lagrangian in question 
contains a massive Yang-Mills part and two higher derivative terms 



C 



AV4> 



-|Tr(F^F L ^ + F R F Rflu ) + m 2 Tr(A^A L " + A R A Rf ") 

- i€Tr(D lt UD v rfF L '" + D ^ D V U F R » V ) + aTnF^ v UF R ^U\ 



where U = exp(2i0/F 7r ), = 0V»/V2, = ±(V^ + = - 

d^A^-dvA^R-ig A^'R, A^ R and the covariant derivative D^U = d^U -igA^U +igU A R . 
The quadratic piece of this Lagrangian is non-diagonal in d^cf) and A^. After diagonalization 
the physical masses are given by 



my 



rrir 



m. 



1 - a 



2 2 

; rn A = m ai 



a 



m + 



4 



(12) 



and F n is related to the physical coupling F n = 135 MeV through 



ZF^ 



1 - 



cr m 



4mg 



(13) 



1 + a m\ 

The couplings <7 aip7r , h aifm and g p7r7r are expressed through g, £, a and the meson masses by 
2Z 2 ( 2 



9 pTTT] 
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2(1 -a) 



i-\(\-z 2 ) 



1 - a 
9i Z 



[l-a)(l-Z 2 ) + 2g& 



l-a) 



Z 2 ) 



(14) 
(15) 
(16) 



Here we retain a non-zero mass for the pion for the purposes of fitting the coupling constants. 
The widths are expressed through these couplings in the following wayQ 



p — >7T7r 



1 |q IV 



and 



ml 



a± — >piT 



127rm?. 



2 / Fp 



I 1 2 i, 
M ft 



a,ip7r 



(17) 



18) 



With the four available parameters g, a, £ and m it is possible to fit both the masses and 
the widths of the p and the a\ [13]. We have refitted these parameters using a recent value 



3 We note that the minus sign in eq. ( Jl8| ) is correct in contrast to refs. jl3], [l4j which are written with an 
incorrect plus sign. 
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of the width, T ai = 400 MeV |15, 16| . There are two possible solutions: 



(A) a = 0.340, £ = 0.446, # = 8.37, 

(B) a = -0.291, £ = 0.0585, # = 7.95. (19) 



which correspond to 



(A) g ai(m = -5.42 GeV, h aifm = -16.7 GeV" 1 , 7 = 0.52 , 

(B) g^pw = 4.25 GeV, h aipn = -2.05 GeV -1 , 7 = 0.33 . (20) 

Here the quantity 7 is the ratio of polar- and axial-vector contributions to radiative pion 
decay. Both solutions are reasonably consistent with the positive experimental value of 



~ 0.4 discussed by Holstein [13]. However it can be shown that the opposite-sign solution, 
(B), is excluded by the QCD sum rule estimates of Ioffe and Smilga ]1^| for the two form 
factors entering the non-diagonal matrix element (ail^Tr). They use couplings g\ and #2 to 
parameterize these form factors in a p-dominance approach, and the relation to g aiP -n and 
is given by 



m 



(21) 



While the absolute values of g\ and g<i obtained in ref. IJT7J contain large uncertainties, they 
are definitely of the same sign, thus ruling out solution (B). Therefore we choose the like-sign 
solution (A). 

We shall display our results for the r.h.s. of eq. (^) multiplied by Q 4 , which according 
to the l.h.s should give the linear relation (2 — b)Q 2 + c. The results from eq. (8) are given 
by the dashed line for the charged pion case in fig. 2. It is seen that there is a good linear 
dependence for Q 2 > 0.9 GeV 2 . We cannot use values of Q 2 larger than plotted in the 
figure since higher states, which are not accounted for, become important and p-dominance 
is not applicable either. There is an excited pion state 7r*(1300) which may contribute for 
the values of Q 2 in question. Its coupling to pre defined through (7T*\p(q)7r(p)) = g* e p -p may 
be roughly estimated using the rather uncertain data [0] on the width, = 200 — 600 
MeV and T n *^. wp = §1^°*. This gives g* ~ 5. The contribution of the it* to the r.h.s. of eq. 
(0) is then 

2g* 2 mi 

- p (22) 

g 2 (Q 2 + m 2 f(Q 2 + ml*) ' { ) 

The result of taking into account the n* is shown in fig. 2 by the full line. It is clear that 
the effect of the ir* is quite small. By fitting a straight line to the curve we estimate b = 1.14 
and c = 1.14 GeV 2 . 
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The matrix element of the gluon field energy density, eq. (||), must be the same for 
charged and neutral pions. This may be used to check our calculation. So, let us now 
consider the case of neutral pions. Isotopic spin invariance forbids the 7r° and a? mesons in 
the s-channel so the lowest allowed state is the u meson. The upn vertex has the the form 

i(Tr(jp)p(q)\u(p + q)) = g^e^ e" a e% Pa q T , (23) 

where and are the polarization vectors of the uj and p mesons. Then the contribution 
of the uj to the r.h.s. of eq. (0) is 

H91 (he)' (24) 

To be consistent we should use the value of the coupling constant g W/97r obtained from the de- 
cay uj — * 7T7 using p-dominance fll8f , g^ pi , ~ 14.9 GeV -1 . The corresponding Q 2 dependence 
of eq. (24) (multiplied by Q 4 ) is shown in fig. 2 by the dashed curve for the neutral case. 

There is, however, an excited state, co>*(1390), which can contribute. The dominant decay 
mode is to the pre channel and taking this to account for the full width of 230 ± 40 MeV 



15], we deduce a coupling constant g^*^ = 5.29 GeV . The result of including both the 



uj and uj* is shown by the full line in fig. 2. There is a noticable curvature and a linear fit 
in this case results in larger uncertainties: 6 = 1 — 1.2 and c = — 0.3 GeV 2 . While the 
value for b agrees with the one obtained from charged pions, it is clear that the intercept 
c is different. This should have been expected since c involves the contributions of quark 
operators and their averages over charged and neutral pions need not be the same. 

Thus, we adopt the value b ~ 1.14, corresponding to a normalization point p ~ Q ~ 1 
GeV. This is in good agreement with the value b = 1.03 deduced from the analysis of ref. 
H in which the matrix element (7r|0"^" d |7r) was extracted from a fit [[0| to the quark and 
gluon distribution functions in the pion. In the leading log approximation the dependence 
on the normalization point is determined by the renormalization group. However, as is well 



known p0[ , operators of the same twist get mixed under renormalization due to radiative 
gluon corrections. The diagonal combinations in the case of two quark flavors are 



Mot 

v 


nu _i nd 

- V + V 


+ e% 


[0] 
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Qu i nd 

- v + V ' 


3nG 






44" 

87 




_ nu nd 

" UV U fIU 








32" 
87 



(25) 

The numbers in parentheses are the anomalous dimensions 7 of the corresponding diagonal 
operators which are renormalized multiplicatively, 

A 7 A «s(/i 2 ) Mg/AgOp) 

Oq = K<O p , K = = - — —— , (26) 

a s {Q 2 ) \og{p Aqcd) 
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where Aqcd ~ 150 MeV. Then the evolution of b, defined by eq. (^), under a change of the 
normalization point is given by 

60i) = ^(l-K 44/mr )+6(g)K 4 ^ 8r . (27) 

It can be seen that according to eq. fl27|) b decreases with // and becomes zero at /i = I.IAqod. 
At the standard normalization point used in QCD sum rules, fi = 0.5 GeV, we get b = 1.06. 
Note that the small value of the normalization point for which 6 = (meaning that there is 
no gluon component in the pion) agrees with the results of ref. ]2l] where it was shown that 
a quark model description of deep inelastic scattering of leptons on nucleons is consistent 
with experimental data provided \i ~ m n . 

Coming back to finite temperatures, the temperature dependence of the condensate (E 2 + 
B 2 ) is determined by the integral over the thermal pion phase space 

where the factor of 3 in front of the integral accounts for the three charged states of pions. 
The structure function Ci in eq. ([!]) is obtained in the same way 



C2(Q,T) = — 2-6+-^ + ... )+0 — . (29) 




where c = \c c harged + \c~neutrai ~ § is the charge averaged value of the constant c. 

Let us now briefly summarize what is known about behavior of condensates at low tem- 
peratures in the chiral limit. The temperature dependence of the usual (Lorentz scalar) 
condensates at low T was considered on the basis of chiral perturbation theory up to three- 
loop order ||. The low T expansion of the quark condensate begins with a term of order 
T 2 /F 2 , because for pions with zero momentum the matrix element (vr|gg|7r) is non-zero and 
proportional to (0|gg|0)/F 2 . In the case of the operator G®G atiV , which is a chiral singlet, 
the one-pion matrix elements vanish. The T dependence of the gluon condensate is related 
through the trace anomaly to (0^)t- The first non-zero contribution to this matrix element 
appears only at the three-loop level. As a result, the T dependence of the gluon condensate 
begins at order T 8 /F£, 

<^l) T = <^> - - l)g (lo^ - i) + . . . , (30) 

where A p ~ 275 MeV is a scale encountered in the three-loop calculation of the pressure of a 
hot pion gas within chiral perturbation theory M. The sign of this contribution corresponds 
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to the melting of the gluon condenate with rising temperature. However, this melting is 
much slower than in the case of the quark condensate, and (G 2 )t is practically constant up 
to T ~ 150 MeV, that is, in the region of applicability of the approximation of a hadronic 
gas. 

One-pion matrix elements of Lorentz non-scalar operators cannot be estimated using the 
soft pion approach, because they are proportional to the pion momentum p. Since p ~ T, the 
corresponding condensates naturally vanish as T — > 0. Since (B 2 — E 2 )t — (B 2 — E 2 )o, we get 
from eq. fl28|) the T dependence of the condensates of chromomagnetic and chromoelectric 
fields, 

(B 2 )t = (B 2 ) + ^T 4 , (E 2 ) T = (E 2 ) + |V , (31) 

where (B 2 ) = — (E 2 ) — 2 x 1CT 2 GeV 4 , using a renormalization scale p = 0.5 GeV. We 
indicate the predicted ratios (B 2 )t/ (B 2 )o and (E 2 )t/ (E 2 )o by the dashed curves in fig. 3. It 
is seen that the T dependence is rather weak at low T and, at T ~ 150 MeV, the condensates 
are changed from their T = value by only about 1%. The fact that the change is small is 
qualitatively consistent with the results extracted from the lattice data ||, however we do 
not agree with the lattice predictions for the sign. We suggest that the lattice calculations 
are probably not sufficiently accurate to predict such small effects. 

We notice that keeping finite would not affect the values of b and c within the ac- 
curacy of our approach. The only differences would appear in the integral over the thermal 
distribution function, eq. (28), and in a lower order contribution to eq. (30). It is straight- 
forward to perform the calculation numerically and this results in the full curves shown in 
fig. 3. We observe that eq. (31) is a good approximation, indeed for the electric field the 
results are indistinguishable. We remark that at very low T, T m n , we have for p = 0.5 
GeV 

(B 2 ) r = (B 2 ) - 0.033m^ 2 T 3 / 2 e- m ^ /T , (E 2 ) T = (E 2 ) + 0.20m^ 2 T 3 / 2 e - m ^ /T (32) 

The numerical effect is exceedingly small, but it is interesting to observe that the magnetic 
condensate (B 2 )^ slightly decreases at very low T before increasing. The behavior of (E 2 )^ 
is, however, monotonic. 

Finally we briefly comment on the effects of higher spin and twist operators. The averages 
of Lorentz non-singlet operators of spin larger than 2 are necessarily proportional to higher 
powers of T and their contribution to thermal correlators will be suppressed by powers 
of T 2 /Q 2 . The operators of spin 2, but of higher twist, are suppressed by p 2 /Q 2 , where 
Ph is some hadronic mass scale ~ Aqcd- In the case of vector currents three twist 4, 
spin 2 operators contribute to the constant c in eq. (|7|) and to disentangle individual 



9 



contributions some extra information must be used. In our opinion, this problem deserves 
further consideration. 

We acknowledge useful discussions with A. Gorski, M. Shifman, E. Shuryak, CS. Song 
and A. Vainshtein. V.E. would like to thank the staff of the Nuclear Theory Group and 
Theoretical Physics Institute, especially Larry McLerran, for the warm hospitality extended 
to him during his stay at the University of Minnesota. This work was supported in part by 
the Department of Energy under contract DE-FG02-87ER40328. 
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Figure Captions: 

• Fig. 1: Diagrams contributing to (a) (7r|0*„ 2 |7r) and (b) (7r|0!f 1 „ 2 |7r). The dashed lines 
correspond to gluons. 

• Fig. 2: The r.h.s. of eq. ([71) multiplied by Q A shown as a function of Q 2 . In the case of 
charged pions, the dashed curve is obtained with tt- and ai-meson intermediate states 
and the full curve also includes the ir* meson. In the case of neutral pions, the dashed 
curve is obtained with an a;-meson intermediate state and the full curve also includes 
the us* meson. 

• Fig. 3: The curves labelled B and E give, repectively, (B 2 ) T /(B 2 ) and (E 2 ) T /(E 2 ) 
as a function of temperature. Normalization point is /i = 0.5 GeV. The dashed curves 
give the results for zero pion mass and the full curves correspond to non-zero pion 
mass. In case E these two curves are indistinguishable. 
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